Mon. Not. R. Astron. Soc. 000, 000-000 (0000) 



Printed 1 October 2008 



(MN WF^ style file v2.2) 



Gamma ray emission from SNR RX J1713. 7-3946 and the 
origin of galactic cosmic rays 

G. Morlino^*, E. Amato^ f, and P. Blasi^'^J, 

QQ ' ^ INAF-Osservatorio Astrofisico di Arcetri, Largo E. Fermi, 5, 50125 Firenze, Italy 

^ Center for Particle Astrophysics, Fermi National Accelerator Laboratory, Batavia, IL, USA 

o ■ 



Accepted . Received 



ABSTRACT 

We calculate the flux of non-tliermal radiations from the supernova remnant (SNR) 
RX J1713. 7-3946 in the context of the non-linear theory of particle acceleration at 
. shocks, which allows us to take into account self-consistently the dynamical reaction 

Qh' of the accelerated particles, the generation of magnetic fields in the shock proximity 

and the dynamical reaction of the magnetic field on the plasma. When the fraction 
\Z^ • of particles which get accelerated is of order ~ 10"'*, we find that the strength of the 

' magnetic field obtained as a result of streaming instability induced by cosmic rays 

. is compatible with the interpretation of the X-ray emitting filaments being produced 

by strong synchrotron losses in lOO^G magnetic fields. The maximum energy of 
accelerated protons is ^ 10^ GeV. If the X-ray filaments are explained in alternative 
" ^ , ways, the constraint on the magnetic field downstream of the shock disappears and 

' the HESS data can be marginally fit with ICS of relativistic electrons off a complex 

Q>^ \ population of photons, tailored to comprise CMB and ambient IR/Optical photons, 

i The fit, typically poor at the highest energies, requires a large density of target photons 

' within the remnant; only a fraction of order ~ 10~® of the background particles gets 

accelerated; the local magnetic field is of order ^ 20/iG and the maximum energy of 
protons is much lower than the knee energy. Current HESS gamma ray observations 
' combined with recent X-ray observations by Suzaku do not allow as yet to draw a 

. definitive conclusion on whether RX J1713. 7-3946 is an efficient cosmic ray accelerator, 

although at the present time a hadronic interpretation of HESS data seems more likely. 
We discuss the implications of our results for the GLAST gamma ray telescope, which 
. should be able to discriminate the two scenarios discussed above, by observing the 

' shape of the gamma ray spectrum at lower energies. 
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1 INTRODUCTION 

The detection of high energy gamma rays from supernova 
remnants has long been considered as the most promising 
tool to test the supernova paradigm for the origin of the 
bulk of cosmic rays observed at the Earth. Several supernova 
remnants (SNRs) have now been detected in gamma rays, 
while u seful upper limits have been obtained for several oth- 
ers (see iTatischeS (|2008l ) and references therein for a recent 
review of observations and their implications). The question 
now to be answered is whether such emissions can be reliably 
considered as a proof of the acceleration of hadronic parti- 
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cles (cosmic rays) or if they should be rather interpreted as 
the result of radiative losses of relativistic electrons. This 
goal can be achieved by investigating the following aspects: 
1) spectrum of the gamma ray emission combined with the 
spectrum of X-rays (and possibly radio emission); 2) high 
energy cutoff of the gamma ray emission; 3) morphology of 
the emission. 

Here we specialize our calculations to the SNR RX 
J1713.7-3946. 

For this remnant, the HESS Cherenkov gamma ray 
telescope has measured a spectrum extending to ~ 100 
TeV, although a power law spectrum with an exponential 
cutoff at g^"^ leads to a b est fit with E^"^ ~ 10 TeV 
lAharonian et al.ll2006l . 120071 '). If this emission is due to de- 
cay of neutral pions produced by nuclear collisions of rela- 
tivistic protons with some target gas in the shock vicinity. 
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then this would imply a cutoff in the spectrum of acceler- 
ated protons at ~ 100 TeV, roughly a factor of 30 short of 
the knee. It is worth keeping in mind that in an expanding 
SNR the highest cosmic ray energy is expected to be reached 
around the beginning of the Sedov phase. Before and after 
that time the maximum energy of accelerated particles does 
not need to be as high as the knee energy. Moreover, differ- 
ent SNRs might accelerate to different maximum energies 
because of the different environment in which they explode. 
In this sense it is only natural that most SNRs do not show 
proton acceleration up to the energy of the knee. 

The strength and topology of the magnetic field in 
the acceleration/radiation region regulate the efficiency of 
particle acceleration (nuclei and electrons) and the radia- 
tion losses of electrons. The rece nt detection of narrow X- 
ray bright rim s in several S NRs (iBam ba et al."2003l, l2005l : 
iLazendic et~al. 2003, 2004; IVink fc L aming 2003) has al- 
lowed to infer the strength of the field, if the thickness of the 
rims is interpret ed as the syn chrotron loss length of the radi- 
ating electrons (lBalletll2006l '). Typical values of 50- 1000^6* 
have been inferred, which suggest efficient magnetic field 
amplification in the shock region. One mor e argument in fa- 
vor of strong magnetic fields was made by lUchivama et akl 
l|2007l ) who interpreted the rapid time variability (of the or- 
ders of few years) observed in some X-ray emitting regions of 
RX J1713. 7-3946 as due to rapid synchrotron cooling. The 
field strength inferred based on this interpretation is as high 
as ~ 1 m G. The same effect was also obs erved in some knots 
of Gas A l|Uchivama fc Aharonianll200§ ). 

There is however a large source of ambiguity in that the 
narrow thickness of the X-ray rims might be due to damping 
of the magnetic field (for instan ce due to wave-wave non- 
linear coupling) (|Pohl et al.ll2005l ) . In this case the magnetic 
field inferred by assuming that the rims are due to severe 
synchrotron losses might be overestimated. We discuss this 
issue at several points throughout the paper. 

Magnetic field amplification could resul t from turbu- 
lent amplification (|Giacalone fc Jokipiil l2007h or from cos- 
mic ray induced streaming instability. An open question is 
whether the large magnetic fields present in the X-ray rims 
are widespread throughout the remnant or are rather con- 
fined in thin filaments with a relatively small filling factor 
(this would be the case if damping plays a role downstream 
of the shock). This difference might have important conse- 
quences on the interpretation of the observed non-thermal 
emission from SNRs, as we discuss below. 

Interestingly enough, the magnetic field strength in- 
ferred from X-ray observations is of the same order of mag- 
nitude as that required by the theory of shock acceleration 
in order to account for particle acceleration up to the knee 
region, provided the diffusion coefficient is Bohm-like. 

From the theoretical point of view, a major improve- 
ment has been achieved in the last few years in that the 
calculations are carried out in the context of the so-called 
non-linear theory of particle acceleration. This theory al- 
lows us to calculate the spectrum and spatial distribution 
of accelerated particles around the shock, their dynamical 
reaction on the shock and the magnetic field amplification 
that cosmic rays induce due to streaming instability. 

Although the non-linear theory has been previously ap- 
plied to SNRs, this paper improves on previous attempts in 
two major aspects: 1) the magnetic field in the shock region 



is calculated from streaming instability and from conserva- 
tion equations instead of being inserted ad hoc in such a way 
to fit the X-ray observations; 2) we include the dynamical 
reaction of the turbulent magnetic field on the shock, which 
is a very impo r tant e ffect in shaping the shock precursor 
(|Caprioli et al.l l|2008al lbl)l. 

The first point is of crucial importance in that it al- 
lows us to achieve a unified picture of the spectra of X-ray 
and gamma ray emission, together with the morphology of 
the X-ray emission, the strength of the magnetic field and 
finally the diffusion properties of the accelerated particles. 
Inclusion of the dynamical reac tion of the field, then, as 
showed bv lGaprioli et al.l (|2008al ). leads to a substantial re- 
duction of the compression in the precursor, thereby also 
reduci ng the concavity of th e spectrum of accelerated par- 
ticles l|CaprioU et al.1 l|2008bD '). A reduced concavity of the 
spectrum should also be expected as a result of the enhanced 
velocity of the scattering c enters when the magne tic field is 
amplified (see for instance IZirakashvili fc PtuskinI ((2000.) ). 

The large values of the magnetic field strength in- 
ferred from X-ray observations and confirmed by our calcula- 
tions clearly favor a hadronic interpretation of the observed 
gamma ray emission, because a subdominant population of 
electrons is sufficient to explain the spectrum and intensity 
of the observed X-ray emission, but is insufficient to produce 
the observed gamma ray emission through inverse Compton 
scattering (IGS) on the GMB and infrared/optical photons 
(IR-l-Opt). The issue of whether the large fields are present 
in large fractions of the remnant volume or are rather con- 
fined in narrow filaments, remains however an open issue. 

For this reason we also investigate the possibility that 
the filaments have some alternative explanation, so that 
there is no constraint on the strength of the magnetic field 
downstream. We still use our non-linear calculations to de- 
rive the spectra of accelerated protons and electrons, but 
force the injection to be low enough to avoid large mag- 
netic fields produced by streaming instability. We find that 
a marginal fit to the gamma ray and X-ray data is possible, 
although a large density of IR photons needs to be assumed 
inside the remnant. The fit to the highest energy part of 
the gamma ray spectrum is not as good as for the hadronic 
interpretation. 

The general conclusion we draw is that although there 
are numerous indications that the observed gamma rays may 
be of hadronic origin, a clear answer on this point may only 
come from: a) the extension of the gamma ray observations 
to lower energies, and b) the detection of high energy neu- 
trinos, whi ch would be produc ed in the same hadronic in- 
teractions (|Morlino et al.ll2008D . 

The paper is organized as follows: in ^ we discuss in 
some detail the general lines of the calculations we carry 
out, including the non-linear theory of particle acceleration, 
the magnetic field amplification and the different channels 
of non-thermal emissions. In !j3] we discuss the results we 
obtain by specializing the calculations to the case of the 
SNR RX J1713. 7-3946. We conclude in g] 

2 DESCRIPTION OF THE CALCULATIONS 

In this section we discuss the technical aspects of the calcu- 
lations we carry out, with special attention for the compu- 
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tation of the spectra of protons and electrons, the magnetic 
field amplification and the radiation processes that we in- 
clude. 



2.1 Spectra of accelerated protons and electrons 

The dynamics of the shock region is dominated by acceler- 
ated protons, while electrons are accelerated in the shock 
environment generated by protons. In this sense the non- 
linearity of the problem is limited to the proton component, 
which also generates the local magnetic field by exciting a 
streaming instability. The spect rum of accelerat e d prot ons is 
computed using the method of lAmato fc BlasH (|2005h . The 
accelerati on time in t he pr esence of a precursor was deter- 
mined bv lBlasi et al.l l|2007l ). and we adopt that calculation. 
It was found there (and later co nfirmed for the parameters 
of RX J1713.7-3946 by lEUisonl (i20_08)) that the precursor 
reduces the maximum energy compared with the naive pre- 
diction of the test-particle theory with the same value of the 
magnetic field upstream of the subshock. 

The maximum momentum is first calculated using as a 
condition the equality between the acceleration time and the 
age of the SNR. However we also illustrate our results in the 
case in which the maximum momentum is determined by the 
finite size of the accelerator. The diffusion coefficient is cho- 
sen to be Bohm-like in the magnetic field generated through 
streaming instability (see §[23J. Since the magnetic field de- 
pends on the location in the precursor, the diffusion coeffi- 
cient is also a function of space: D(p,x) — {l/'i)p(? / {eB{x)). 
The two conditions mentioned above read: 



~ tsNR , 



Uo 



= aRsf 



(1) 



where a < 1 is the fraction of the shell radius (Rsnr) where 
the escape of particles at p ~ pp , max occurs. 

Although the basic st ructure of the calcu lation 
is the same p ropos ed by lAmato fc BlasH l|2005l ) and 
lAmato fc BlasH (|2006h . the crucial new aspect taken into 
account here is the dynamical reaction of the self-generated 
magneti c field. We introduce t his e ffect following the treat - 
ment of ICaprioh et al] (j2008al ) and ICaprioh et al.l (|2008bh : 
the conservation equations at the shock and in the precur- 
sor are modified so as to include the magnetic contribution. 
The compression factor at the subshock, Raut, and the total 
compression factor, Rtot, are deeply affected by this change, 
in that the ratio Rtot/Rsub (the compression factor in the 
precursor) , decreases when the amplified magnetic field con- 
tributes a pressure which is comparable with the pressure of 
the thermal gas upstream. In turn this reflects in spectra of 
accelerated particles which ar e closer to power laws , though 
a concavity remains visible (|Caprioh et all l2008bl ). as we 
discuss below. 

The normalization of the proton spectrum is an output 
of our non linear calculation, once a recipe for injection has 
been established. Following Blasi et al. (2 005 ). particles are 
injected immediately downstream of the subshock. The frac- 
tion of gas particles crossing the shock from downstream to 
upstream, jjinj, can be written as 

»7inj=4/(3y¥)(i?,,b-l)C3e-«'. (2) 

Here ^ ~ 2 — 4 is defined by the relation pinj — 5pth,2, where 
Pth,2 is the momentum of the thermal particles downstream. 



5 parametrizes the poorly known microphysics of the injec- 
tion process, but one should stress that pth,2 is an output of 
the problem: as a result, the injection efficiency is affected by 
the dynamical reaction exerted by the accelerated particles 
and by the amplified magnetic field. 

The spectrum of accelerated electrons at the shock, 
fefiip), is easy to calculate for p < Pe.max, with pe,max the 
maximum electron momentum. In fact the slope of the elec- 
tron spectrum at given momentum p is the same as the slope 
of the proton spectrum at the same momentum, assuming 
that both electrons and protons experience the same diffu- 
sion coefficient. The relative normalization of the two spec- 
tra is unconstrained a priori, while the normalization of the 
proton spectrum, as stressed above, is an output of our non 
linear calculation. 

The maximum energy of electrons is determined by 
equating the acceleration time with the minimum between 
the time for energy losses and the age of the remnant. The 
loss time of electrons over a cycle of shock crossing needs 
to be weighed by the residence times upstream and down- 
stream, therefore the condition for the maximum momen- 
tum, in the loss dominated case, can be written as: 



tacc{p) = [Ti,l(_Bl,p) + Tl,2{B2,p)] 



tr,l{p) 



+ 



U.2{p) 



n.i{Bi,p) Ti, 2(52, p) 



where r; denotes the loss time, and the indexes "1" and 
"2" refer to quantities measured upstream and downstream 
respectively. The residence times in the context of non lin- 
ear theory of particle accel eration can be writ ten explicitely 
(from Eqs. (25) and (26) of iBlasi et all (|2007l ')'). Eq. ^ can- 
not be solved for Pe,max analytically, contrary to the case of 
acceleration in the test particle regime. The numerical solu- 
tion is however easy to obtain with standard techniques. 

When only synchroton losses are important, we can pro- 
vide an approximate expression for the solution of Eq. Q. 
We use the function Up{p) defined as the mean plasma ve- 
locity that a particle with mo mentum p ex periences in the 
precursor region (see Eq. (8) in lAmato fc Blas i (2005)). The 
function Up{p) = Up{p)/uo is such that Up{p) < 1 for any 
momentum. The electron maximum momentum can then be 
written in an implicit way as: 



Pe,max 



mic* uo 



2 V eSiro 



Up (Pe.max) 



^ ^tot Up (Pe.max) 
1 -I- RnRtotUpiPe , max J 



where ro is the classical electron ra dius. One may noti c e that 
Eq. (gl) reduces to the Eq. (11) of iBerezhko fc Volkl (|2006t ) 
if one poses Up — 1. Eq. (|4]) leads to an error of less than 
10% with respect to Eq. (jSJ (but we stress that the results 
presented in this work are all obtained using Eq. ((3])). 

The spectrum of electrons at energies around and above 
Pe.max, namely the shape of the cutoff is harder to calculate 
in the context of a fully non linear calculation. In the test 
particle case and f or strong shocks th e solution has been cal- 
culated bv lZirakashvili fc AharonianI |2007 '). Since the spec- 
tra we find for electrons at p < Pe,max are not far from being 
power laws w ith slope ~ 4, we adopt the modi fication factor 
calculated bv lZirakashvili fc AharonianI l|2007l ') and we write 
the electron spectrum at the shock in the loss dominated 



fe,0{p) = Kep fp,0{p) 1 + 0.523 (p/pe,max)'l "e"''"/''^ 



.(5) 
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The most important aspect of this expression is the fact that 
the cutoff is not a simple exponential, which in turn reflects 
in the shape of the secondary emission radiated by electrons. 
This important conclusion differs from the numerous calcu- 
lations previously carried out by other authors, where an 
exponential cutoff was assumed for the electron spectrum. 
On the other hand, if the maximum momentum is indeed de- 
termined by the age of the remnant, then the cutoff shape is 
expected to be exponential. The constant Kep accounts for 
the different normalization between the electron and proton 
spectrum. 



2.2 Magnetic field amplification and compression 

The turbulent magnetic field close to the shock can be en- 
hanced by several physical processes. Here we concentrate 
on the amplification due to streaming instability induced by 
cosmic rays accelerated at the shock and we show that this 
mechanism allows us to reproduce the strength of the mag- 
netic fields which are inferred by assuming that the thickness 
of the X-ray filaments is due to severe synchrotron energy 
losses. On the other hand, if the filaments are interpreted 
in a different way, very low efficiencies of proton accelera- 
tion are required in order to avoid effective magnetic field 
amplification. 

The streaming of cosmi c rays upstream of a shoc k ex- 
cites both resonant (iSkillingl (| 19751 )) and non-resonant (|Belll 
|2004) modes. The latter modes are expected to dominate 
for fast shocks, efficiently accelerating cosmic rays, most 
likely during the free expansion phase of the post-supernova 
evolution and at the very beginning of the Sedov phase 
l|Amato fc Blasi 2008). At any later time, resonantly excited 
Alfven waves should account for most of the field ampli- 
fication in the shock vicinity. Although the predictions of 
quasi-linear theory should be taken with caution when con- 
sidering the regime of strongly non-linear field amplification, 
SB/B ^ 1, at present no better treatment is possible. The 
strength of the magnetic field at the position x upstream, 
6B(x), in the absence of damping, can be estimated from 
the saturation conditi on, that, for modified shocks, reads 
l|Caprioh et al.l (|2008bl )): 



Pw{x) = U{x) 



-3/2 



Pw,0 + 



U{x) 



4. Mam 



(6) 



where Pm{x) — SB{x)^ /{STrpoUo) is the magnetic pressure 
normalized to the incoming momentum flux at upstream in- 
finity. In Eq. (|6l), U{x) = u{x)/u{o) and Ma,o = uq/va is 
the Alfvenic Mach number at upstream infinity where only 
the background magnetic field. Bo, assumed to be parallel 
to the shock normal, is present, p^^ accounts for the mag- 
netic pressure due to a pre-existing magnetic turbulence. 
This term is usually set to zero in the following calculations, 
unless otherwise specified. 

Eq. ((Ojl correctly describes into account the effect 
of compression in the shock precursor through the term 
U{x)~^^'^. For the upstream temperature that we adopt in 
RX J1713. 7-3946 (see below), damping in the upstream re- 
gion is expected to be negligible. 

This expression has been derived for Alfven waves, 
therefore the perturbation is perpendicular to Bo, namely 



in the plane of the shock. The magnetic field downstream of 
the subshock is further enhanced by compression: 



B2 — Rs u b Bi , 



(7) 



where Bi is the magnetic field immediately upstream of the 
subshock. 

The combined action of cosmic ray induced magnetic 
field amplification and field compression at the subshock re- 
sults in a situation in which the magnetic pressure may easily 
exceed the pressure of the background plasma. As shown by 
(|Caprioh et al.ll2008ah . this leads to a reduction of the com- 
pressibility of the plasma, so that the value of Rtot drops 
compared with the cases in which this dynam ical reaction of 
the magnetic field is ignored. As stressed bv lCaprioli et af] 
(2008i) this effect alleviates the need for invoking poorly 
known processes such as the so-called turbulent heating. 



2.3 Radiation produced by accelerated particles 

In this section we briefly describe the calculations of the 
fluxes of nonthermal radiation produced by electrons and 
protons in the SNR. We concentrate on the following chan- 
nels: 1) gamma ray production by generation and decay of 
neutral pions in inelastic pp collisions; 2) synchrotron emis- 
sion of electrons; 3) ICS of electrons on the CMB and in- 
frared/optical photons in the remnant. 

The flux and spectrum of gamma rays from tt" decay are 
cal culated following th e analytical approximations discussed 
bv lKelner etU] l|2006l ). 

The cutoff in the gamma ray spectrum resulting from 
-K^' decay can be written as 



^y.cut ^ 0.15 -Ep.max , 



(8) 



where 



is the maximum energy of accelerated pro- 



tons. Qualitatively, the spectrum of gamma rays at energy 
E <C i?^,cut reproduces the shape of the parent proton spec- 
trum. The gamma ray spectrum is therefore expected to 
show some level of concavity, typical of acceleration at cos- 
mic ray modifled shocks. At energies lower than the pion 
mass, the gamma ray spectrum is found to rapidly drop. 

Synchrotron emission is produced by accelerated elec- 
trons in the magnetic fleld amplified by streaming cos- 
mic ray protons. The calculation of the synchrotron emis- 
sion is carried out by usin g the exact synchrotron kernel 
dRvbicki fc LightmanI [l985l ) . The synchrotron spectrum is 
cut off at a photon energy 



E, 



syn,cut 



3 X lO^'^Bi 



/ Ee,7nax \ 

\ TeV J 



keV, 



(9) 



where Bioo is the local magnetic field in units of lOOfiG. It 
is easy to check that for Ee,max ~ 10 TeV, the synchrotron 
emission is cut off in the keV energy range. It is also worth 
keeping in mind that in the straightforward case of Bohm 
diffusion at a linear shock, the maximum energy of elec- 
trons scales as B~^^^, therefore Esyn,cut is independent of 
the magnetic field. In our calculations, the non-linear effects 
induce a weak dependence on magnetic field, as discussed in 

§13 

In this paper we compute the synchrotron radiation 
spectrum limited to X-ray band. The radio emission, be- 
ing produced by lower energy electrons, with much longer 
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lifetimes, is affected by the temporal evolution of the super- 
nova remnant, which strongly depends on the environment 
in which the blast wave propagates. 

The X-ray emission is confined to a narrow region be- 
hind the shock, whose thickness is determined, depend- 
ing on the radiating electron energy, by either the diffu- 
sion or advection length covered by electrons during their 
loss time. As we show at the end of next section the con- 
tribution to the X-ray radiation from the upstream re- 
gion is negligible, due to the much lower magnetic field 
(see Fig. [5] and the corresponding discussion) . The thick- 
ness of the radiating rim can be estimated by solving the 
diffusion-convection equation in the downstream plasma: 
U2{dfe/dx) — Dd^ fe/dx^ — /e/Tsyn, where B2 has been as- 
sumed to be spatially constant (i.e. we are neglecting any 
magnetic damping mechanisms). The solution has the form 
f{p,x) (X exp(— x/Ai?rini(p)) whcrc Ai?rim(p) is the spatial 
scale of the emit ting region for electron s with momentum p 
and is given by (|Berezhko fc Volkll20o4 ): 



2D2/U2 



- f 



(10) 



The observed thickness of the rim is obtained by correcting 
the previous expression for projection effects. In the ideal 
case of a spherical shock, an exponential pro file translates 
into a projected thickness AJJobs — 4.6AJ?rim (jParizot et al.l 
I2OO6I ). where A7?obs is defined as the size of the region over 
which the observed brightness drops to half of the maximum. 

Finally, electrons contribute to gamma ray emission 
through their ICS of ambient photons. The main contribu- 
tion to gamma ray emission comes from upscattered dust- 
processed infrared photons , having a blackbod y spectrum 
with temperature ~ 20 K (jSchleeel et al.lll998h . Following 
the instance of the IR-f Opt photon background in the inter- 
stellar medium (ISM) we assume that the ratio of the optical 
to infrared energy densities remains ~ 20, while the energy 
density of IR light, Wir, is left as a free parameter (in the 
ISM, Wir ^ 0.05 eVcm-^). 

The flux and spectrum of ICS photons is calculated by 
using the exact ker nel for ICS, with the full K lein-Nishina 
(KN) cross section (jRvbicki fc LightmanlllOSSh . The spec- 
trum is cut off at an energy 



Eics.cut 



\ Tev J V^y 



TeV 



(11) 



if Ee,niax ^ 260(eph /eV) " ^ GeV and Eics,cut = -Ec!,max oth- 
erwise. The boundary between the two regimes corresponds 
to the setting in of the KN regime. For the photon back- 
ground we adopted here (eph « 2 x 10~^ eV), the threshold 
is at energy ~ 150 TeV, therefore for the case of RX J1713.7- 
3946, the cutoff region falls well within the Thompson regime 
and Eq. applies. 



3 RESULTS: RX J1713. 7-3946 

In this section we describe the results of our calculations 
for the case of SNR RX J1713.7-3946. A major difficulty 
in carrying out detailed calculations on this remnant is 
that neither its dist ance nor its age are well established. 
iKovama et al.l (|l997l ) infer a di stance of ^ 1 kpc, based on 
ASCA X-ray observations, while [Slane et all (|l999l ) estimate 



a distance of ~ 6 kpc. Here we adopt the value of 1 kpc, 
which appears to be supported by some more recent work 
(|Fukui et al.ll2003l : iMoriguchi et aLlbOOSl ) based on the spa- 
tial correlation between the remnant and a molecular cloud 
located at that dist ance. This estimate was also indepen- 
dently confirmed bv ICassam- Chenai et al] l|2004D based on 
the comparison between the interstellar absorption as mea- 
sured using X-ray and radio observations. Since the observed 
angular size of the remnant is ~ 60', the remnant radius can 
be estimated to be -Rsnr — 10 pc, which would lead to es- 
timate an age of ~ 1600 yr, consistent with t he historical 
Chinese record of a SN exploded in AD 393 (|Wang et al.l 
Il997f ). 

Several authors, based on different arguments, have 
reached the conclusion that RX J1713. 7-3946 is a core col- 
lapse SN of type Il/Ib, expanding into a hot and very 
diluted bubble created by the wind of a massive progen- 
itor star with initial mass between 15 < M < 20 Mq 
dBerezhko fc Volkll2006l:ISlane et al.|[l999l : iFukui et al.ll2003l : 
ICassam-Chenai' et al.l 20041 ). These authors also conclude 
that the bubble must have a radius of 15 ^ Rb 20 pc, 
hence the shock front of the remnant is most likely still prop- 
agating into the low density shocked wind, except perhaps 
for the SW region where the shock seems to impact the 
cold and dense shell of unshocked material. This scenario 
lead s to a tem perature for the bubble typically jj 10^ K 
( Ch evalier fc Liang lOSg.). In order to account for the inho- 
mogeneity of the medium, we adopt here a somewhat lower 
value of To = 10^ K. We checked however that our results 
are very weakly dependent on the value of the gas tempera- 
ture at upstream infinity, as long as it is larger than ~ 10* 
K. 

Other uncertainties that affect the results of our calcula- 
tions concern the density of the background medium, no , and 
the strength of the background magnetic field (before ampli- 
fication takes place). Here we adopt no = 0.12 cm"'', which 
is somewhat dictated by the fitting to the HESS data. Some- 
what larger values of no could be still compatible with HESS 
observations if a fudge factor is introduced to account for the 
fact that the gamma ray emission is not homogeneous over 
the SN shell. Several attempts at measuring the ambient 
gas density have been made. An upper limit of ~ 2.6 cm"'' 
comes from the molecular hydrogen column density in parts 
of the SNR without detect able CO emi ssion that can be de- 
rived ba sed on NANTEN (iFukui et al.ll200a ') measurements 
(see also I Aharonian et al.l 1 20061 ) for a discussion) . 

The other existing constraint was inferred by 
ICassam-Chenai' et al.l l|2004l ) based on the lack of thermal X- 
ray emission in XMM data. The limit reads no < 0.02 cm"''. 
However, it is worth keeping in mind that this limit was de- 
rived by assuming that electrons and protons are in thermal 
equilibrium downstream, but the same authors suggest that 
this is likely not the case. They also investigate the non- 
equilibrium case but apparently they used an equal temper- 
ature of electrons and protons in this second case as well. 

An upper limit of 4500 Km/s to t he ve locity of the 
ejecta was inferred by lUchivama et al. I l|2007h based on a 
direct measurement of the transient X-ray emitting region. 
On the other hand a value uq < 1500 Km/s would not allow 
for acceleration of electrons up to high enough energy to 
explain the observed X-ray emission. Here we adopt a value 
of uo = 4300 Km/s necessary to fit the cutoff in the X-ray 
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Figure 1. Cutoff energies for photons produced by synclirotron 
emission {dotted line, in keV), pion decay (solid) and IC scatter- 
ing {dot-dashed for scattering on IR radiation, dashed for CMB), 
as a function of the downstream magnetic field B2 ■ Also shown on 
the upper x-axis of both panels is the corresponding value of the 
background magnetic field Bq that was used for the calculation. 
Notice that while the B2-scale is logarithmic, the i?o-scale is arbi- 
trary and built to show the correspondence between Bq and B2. 
The upper panel shows the case where Bq is resonantly amplified 
by the streaming instability. Here the thick vertical lines indi- 
cate the solutions having the X-ray rim thickness equal to 1%, 
2% and 3% of the SNR radius respectively. In the lower panel 
we assume no magnetic field amplification, but simple compres- 
sion of Bo in the precursor and at the shock. The results in both 
panels are obtained with ^ = 3.7 and no = 0.12cm~'^, and with 
the maximum proton momentum determined by the condition 

tacc{pmax) = 1600 yr. 



spectrum as detected by Suzaku. A lower value of about~ 
3000 Km/s would be needed to fit the ASCA data, but we 
do not discuss tliis case any furtfier here, given the superior 
quality of the Suzaku data. 

The model parameters are further constrained by the 
fact that the magnetic field in the downstream plasma, B2, 
is an output of our calculations, for given So, as it is the 
result of cosmic ray induced magnetic field amplification 
and further enhancement by compression in the precursor 
and at the subshock. The correct solution of the problem 
is found by not only requiring that the fiuxes of observed 
radiations are fit, but also requiring that B2 corresponds to 
the value inferred from the measured thickness of the X-ray 
rii n. For SNR RX J1713 . 7-3946, we use the data presented 
in iLazendic et al.l (|2004h where the authors use a CHAN- 
DRA observation at 5 keV to infer the rim thickness in the 
northwestern region of the SNR. The procedure leads to 
Ai?obs ^ 0.02i?sNR (see their Fig. 9). Clearly this thick- 
ness does not lead to an estimate of B2 if due to dam ping of 
the magnetic field (as proposed bv lPohl et al.l l|2005l )') rather 
than to particles' energy losses. Below, we first assume that 



the downstream magnetic field can indeed be derived from 
the extent of the X-ray filaments and then we investigate 
the consequences of relaxing this constraint. 

A first indication of the origin of the observed radiations 
can be obtained by looking at the cutoff frequencies in the 
different bands. In the upper panel of Fig. [1] we show the 
cutoff energies for synchrotron emission (dotted line, in units 
of keV), tt" decay (solid line, in TeV) and ICS (dash-dotted 
line for ICS on the infrared light and dashed line for ICS 
on the CMB) as given by Eqs. (O, (HJ and Hi}. The plot 
is obtained for fixed ^ = 3.7 and varying Bo between 0.1 
and 10 nG (upper a:;-axis) , while the lower x-axis shows the 
corresponding value of B2- The vertical thick solid line shows 
the solution that provides the central value of the measured 
rim thickness (corresponding to B2 ~ 100 fiG) while the 
dashed vertical lines bound the allowed region corresponding 

to 1% < Ai?obB/7?SNR < 3%. 

For the range of values of Bo considered here, it is ap- 
parent that ICS cannot provide a high enough cutoff energy 
so as to explain HESS data. The latter require £7, cut ~ 10 
TeV, while the highest value of the ICS cutoff energy ob- 
tained for scattering on the IR light in the presence of an 
extremely low background magnetic field (Bo ~ O.l^G) is 
~ 4 TeV. The corresponding cutoff energies for ICS off the 
CMB photons are lower by roughly the ratio of tempera- 
tures: ~ Tdust/TcMB ~ 7. The general trend of the ICS 
cutoff energy is to decrease as a function of Bo (and B2) 
because larger values of the magnetic field lead to stronger 
synchrotron losses and thereby lower maximum energies of 
the electrons. There are however two subtleties. First, taking 
into account nonlinear effects in the computation of the ac- 
celeration time reduces somewhat the dependence of Ej^cut 
on the magnetic field, with respect to the expectations of the 
linear shock analysis. Moreover, for values of B2 lower than 
~ 20 ^G, the maximum electron energy is determined by the 
finite age of the remnant rather than by synchrotron losses. 
Hence, for values of the magnetic field below this thresh- 
old, the maximum electron energy increases with increasing 
magnetic field (because of the more effective diffusion) and 
reaches a maximum for B2 ~ 20 fiG. This corresponds to 
the highest possible value of E-y^cut- 

It is important to stress that the conclusion that the 
cutoff frequency for ICS falls short of explaining HESS data 
by about one order of magnitude is independent of argu- 
ments related to the normalization of the flux, which we 
comment upon below. 

The cutoff energy for gamma rays produced by tt" decay 
is > 10 TeV provided Bq > (0.3— 1) ^G, which interestingly 
enough is also the region where B2 is consistent with the 
thickness of the X-ray rims, if these are interpreted as the 
result of severe synchrotron losses downstream. 

The cutoff energy of synchrotron emitted photons, 
which, assuming Bohm diffusion, in linear theory would not 
depend on the value of the magnetic field, now acquires a 
weak dependence on Bo, and varies between ~ 1 and 2 keV 
for the considered range of magnetic field values, granted 
that B2 ~ 20 fiG. Things change drastically when B2 is 
lower than this threshold and the maximum electron en- 
ergy is determined by the competition between acceleration 
time and age of the system, rather than synchrotron losses 
(lower panel of Fig. [1} . 

The adopted value ^ = 3.7 corresponds to a fraction 
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Figure 2. The thick lines in the lower panel show the cutoff 
frequencies for the radiative processes (just as in Fig. [TJ now 
computed as a function of the injection efficiency for fixed Bq = 
2 iiG. The thin soUd fines show the cutoff frequency for tt^ decay 
when the maximum proton energy is computed according to the 
definition on the right of Eq. ((TJ with a = 0.08 and 0.15. The 
downstream magnetic field, resulting from the amplification and 
further compression of Bq, is shown in the upper panel, where 
the vertical lines have the same meaning as in Fig. [T] The upper 
X-axis of this panel displays the fraction of injected particles in 
units of 10~^. 



of injected particles jyinj = 1.2 • 10~*. The resulting acceler- 
ation efficiency, e, defined as the percentage of the energy 
flux crossing the shock that ends up in relativistic parti- 
cles, varies from ~ 68% (for Bq =0.1 ^G) to ~ 36% (for 
Bo = 10 /iG). 

It is very important to notice that the solutions cor- 
responding to the curves in the upper panel of Fig. [T] are 
characterized by moderately low values of the compression 
factor in the precursor {Rtot ~ 12), despite the high accelera- 
tion efficiencies. This result is purely the consequence of the 
dynamical reaction of the amplified magnetic field on the 
shock. No turbulent heating has been introduced in these 
calculations. One may wonder how the results would change 
if no cosmic ray induced magnetic field amplification were 
present, so that the magnetic field in the precursor and at 
the shock would only be produced by compression of a tur- 
bulent field at upstream infinity, with SB ^ Bo- In this case 
the turbulent field (assumed to be perpendicular to Bq) is 
compressed according to B±{x) = Bq{uo/u{x)). In the bot- 
tom panel in Fig. [l] we show the cutoff energies obtained in 
this situation (we still assume ^ — 3.7). For 0.1 < Bo < 10 
liG, the total compression factor ranges between 14 and 7, 
and B2 reaches up to ~ 100 fiG. Remarkably, the highest 
value of the maximum energy of gamma rays produced by 
ICS (off the IR photons) is ~ 7 TeV and is reached for 
Bo ~ 2fiG {B2 ~ 20fiG). Lower and higher values of Bq 
lead to smaller cutoff energy for the ICS emission. 

The results of the calculations are very sensitive to the 
recipe adopted for the injection of cosmic rays from the ther- 
mal pool. As discussed in § 12.11 we model the inj ection with 
the th ermal leakage-like model introduced by iBlasi et al.l 
(|2005l '). In Fig. [2] we illustrate the dependence of B2 (up- 
per panel) and of the cutoff energies (lower panel) on the 
injection parameter ^, varying between 3.4 and 4.1 (corre- 
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sponding to 6 x 10"* > ryinj > 8 x 10""). The upstream 
magnetic field is _Bo = 2 /iG. The curves are labelled as in 
Fig. [T] with the addition of two thin solid lines representing 
the cutoff energy for tt" decay when the maximum energy of 
the parent protons is calculated with a = 8% and a = 15% 
in Eq. ((T}. The vertical thick solid line in the upper panel 
shows the solution providing the central value of the mea- 
sured rim thickness (corresponding to B2 = 100 fiG) while 
the dashed vertical lines bound the allowed region corre- 
sponding to 1% < Ai?obs/-RsNR < 3%. 

The cutoff energy for gamma rays due to ICS is again 
lower than ~ 10 TeV, with the possible exception of ICS on 
IR light for the highest values of ^. For instance, for ^ ~ 3.9 
the cutoff energy for ICS is ~ 12 TeV, but in this case the 
cutoff frequency for pion decay is at the same level. In this 
limit case the issue becomes the absolute normalization. In 
any case such situation also corresponds to B2 ~ 50/iG, 
quite lower than inferred from X-ray measurements (com- 
pare with the region between the two vertical dashed lines 
in the upper plot). 

The above discussion leads to conclude that there is 
some degeneracy in the choice of the values of the param- 
eters of this problem. If we try to fit the spectral energy 
distribution in the X-ray and gamma ray bands and the 
thickness of the X-ray filaments simultaneously, we find 
that the best fit is obtained for 1 < Bo < 3 fiG and 
3.7 < ^ < 3.85. The resulting acceleration efficiency varies 
in the range 20% < e < 60%. In Fig. [3] we show a possible 
solution with Bo = 2.6/iG and ^ = 3.8. The contribution of 
tt" decay to the gamma ray flux is plotted as a thick solid 
line. The solid curve in the lower energy part represents the 
synchrotron emission of electr ons, superimposed on the re- 
cently published Suzaku data l|Tanaka et al.ll2008l ) . The flux 
of ICS photons when the target is made of CMB photons 
is plotted as a dashed line. The ICS on the IR-f Opt pho- 
tons, assumed to have a density equal to the average value in 
the Galaxy (1 eV cm~^) corresponds to the lower thick dash- 
dotted line. For illustration purposes, we also plot the contri- 
bution of ICS on a IR-|-Opt background with energy density 
200 times larger (higher dash-dotted line). It is worth notic- 
ing that while synchotron emission is always the dominant 
loss mechanism for electrons, the rate of energy losses due 
to inverse Compton scattering mainly comes from interac- 
tion with the optical photons. On the contrary, the leptonic 
gamma ray flux in the HESS band is mostly contributed by 
upscattering of the IR photons. 

Fig. [3] shows a few important points: 1) in terms of 
shape, normalization and position of the cutoff the contri- 
bution of tt" decay provides a satisfactory flt to HESS ob- 
servations; 2) the ICS on the CMB photons is characterized 
by too low a normalization and too low cutoff energy; 3) for 
ICS on the IR-|-Opt light, the normalization can be raised 
to the observed flux level only if a very large photon flux is 
assumed; 4) even in this case the position of the cutoff in 
the gamma ray spectrum turns out to be below the observed 
one; 5) while being below the sensitivity limit of EGRET, 
the predicted fluxes should be observable by GLAST. 

As the HESS observations show l|Aharonian et all 

[20o3), the gamma ray spectrum in the TeV region can be flt 
by a cutoff with shape exp[— (i?/Bo)^ ] with /3 close to 0.5. 
As shown here (see also iKelner et all (|2006l )). this result is 
perfectly compatible with the spectrum of gamma rays from 
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Figure 3. Spatially integrated spectral energy distribution of RX 
J1713. 7-3946 in the hadronic scenario, for no = 0.12 cm~'^, uq = 
4300 km/s, Bq = 2.6 ^G, ^ = 3.8. The following components are 
plotted: synchrotron (dotted line) thermal emission (dotted line), 
Compton scattering with CMB (dashed line) and with Opt+IR 
background (dot-dashed line). The contribution from pion decay 
is shown as a thick solid line and corresponds to pp,max = 1.26 X 
10^. HESS data taken from 2003 to 2005 are plotted together 
with Suzaku data in the X-ray band. Also EGRET upper limit 
and GLAST sensitivity are shown. 
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Figure 4. Spatially integrated spectral energy distribution of 
RX J1713. 7-3946 in the leptonic scenario. The following compo- 
nents are plotted: synchrotron (thin solid line) and thermal (thick 
dotted line) electron emission, ICS component for CMB (dashed 
line), optical (dotted line) and IR (dot-dashed line) photons, and 
the sum of the three (thick solid line). The Opt-|-IR components 
are assumed to have energy density 24 times the mean ISM value. 
The other parameters are: no = 0.01 cm"'', uq = 4300 km/s. 
So = 1.5 fiG and ^ = 4.1. The experimental data are the same as 
in Fig. [3] 



the decay of neutral pions, if the cutoff in the spectrum of 
the parent protons is exponential. 

A few words are in order concerning the maximum mo- 
mentum of protons. The value of Pp,max that provides the 
best fit to the data is 1.26 x 10^ m^c, which is obtained when 
the maximum proton energy is calculated by using the con- 
dition of finite size of the accelerator with a — 0.08 in Eq. ([1]) 
(see also Fig. [2]). The condition on the age of the remnant 
would have resulted in a maximum proton energy larger by 
a factor ~ 2. 

A criticism to the scenario just discussed could be raised 
in that the fiux of thermal X-rays plotted in Fig. [3] as a thin 
dotted line exceeds the Suzaku observations. This is how- 
ever true only if protons and electrons in the downstream 
plasma are in thermal equilibrium, w hile there are reasons 
to bel ieve that this is not the case (jCassam-Chenai et al.l 
l|2004h l. If one assumes Te = 0.01 Tp the absence of thermal 
X-ray emission is naturally explained (thick dotted line in 
the same plot). 

The Suzaku data points are fit extremely well by 
our calculations: a fact which at least confirms that the 
shape of the electrons spectrum close to the cutoff is cx 
exp[—{E / Ee,max)^] rather than a simple exponential. In 
other words, the maximum energy of electrons in this rem- 
nant is determined by energy losses, rather than by the finite 
age of the system. Although not very stringent, this argu- 
ment suggests that the magnetic field cannot be too low, or 
otherwise ICS losses would have to be dominant, implying 
the presence in the remnant of a density of IR photons much 
larger than in the ISM. 

What if the narrow X-ray filaments are due to mag- 
netic field damping or other processes other than severe 
synchrotron losses? In this case there is no constraint on 
the strength of the magnetic field at the shock, therefore 
one may try to test the possibility of fitting the gamma ray 
data with the ICS emission of high energy electrons. We 
carry out the previous calculations, including the magnetic 
field amplification due to streaming instability but choosing 
now a large value of the parameter ^, so as to reduce the in- 



jection. This leads to two main effects: 1) lower efficiency of 
proton acceleration and 2) correspondingly lower efficiency 
of magnetic field production. 

For the same values for the supernova as above (espe- 
cially the distance and velocity of the shock), and assuming 
no = 0.01 cm~'^ and ^ = 4.1, we obtain the results illus- 
trated in Fig. 2] These parameters correspond to a fraction 
of accelerated particles of order ~ 7 x 10~^ and a cosmic ray 
energy conversion fraction of ~ 2%. The curves are labelled 
as follows: the solid and dotted line on the left side of the 
plot are the syncrotron and thermal electron emission re- 
spectively. The thin dashed, dotted and dot-dashed curves 
are the ICS on the CMB, optical and IR photons respec- 
tively, while the solid thick line is the sum of the three. The 
assumed energy density for the optical and IR radiation is 
24 and 1.2eV/cm^ respectively. 

From Fig.|4]we conclude that a fit to HESS data on RX 
J1713. 7-3946 can be obtained within an ICS scenario at the 
price of: 1) assuming an IR background in the SNR which 
exceeds the ISM value by a factor ~ 24; 2) neglecting the 
highest energy data points of HESS. 

The first issue has to do with the plateau-like shape 
of the gamma ray spectrum in the lower energy part, which 
cannot be fit unless a second component of the photon back- 
ground (at highe r frequencies) i s assu med. This argument 
was also made bv lTanaka et al.l (|2008l ). 

The second point is due to the fact that despite the 
small efficiency for proton acceleration, the magnetic field in- 
ferred for the downstream plasma is still ~ 20/iG, therefore 
the maximum energy of electrons is still not large enough 
to have gamma ray emission above 10 TeV. Moreover it 
is worth noting that, in spite of the high energy density 
assumed for the optical background, i.e. 24 eV/cm~^, the 
synchrotron losses still dominates over the IC losses because 
the scattering occurs in the KN limit. As a consequence the 
spectrum of electrons is cut off as exp[—{E / Ee,max)^], re- 
fiecting in a gamma ray spectrum with a rather sharp cutoff 
which does not extend to sufficiently high energies to reach 
HESS highest energy data. This is a point which is often 
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Figure 5. Integrated spectral energy distribution for the leptonic 
scenario as in Fig. |4] Here only the synchrotron emission and 
the ICS for CMB are shown, separated in the contributions from 
the upstream (dot-dashed line) and from the downstream region 
(dashed line). The solid line is the total contribution. 



overlooked in papers investigating the ICS interpretation of 
the HESS data. 

The fluxes shown in Figs. [3] and H] are obtained by inte- 
grating over the emission volume, comprising the upstream 
and downstream regions. It is instructive to show the emis- 
sion in the different bands from the upstream and down- 
stream regions separately, since it illustrates some of the non 
linear effects discussed above. In order to do that, we concen- 
trate on the case of inefficient acceleration — 4.1) which 
leads to the leptonic scenario. We stress again that even in 
this case the magnetic field is generated upstream by the 
streaming cosmic ray protons and further compressed in the 
precursor and at the subshock. Moreover, for simplicity we 
use only the CMB photons as a target for ICS (this is the rea- 
son why the predicted fiuxes do not fit HESS data points). 
The fiuxes from downstream (dashed lines) and upstream 
(dot-dashed lines) are illustrated in Fig. \S\ The solid lines 
represent the total emission, as the sum of the upstream and 
downstream contributions. For both synchrotron emission 
and ICS the main contribution comes from the downstream 
region, at almost all photon energies. For synchrotron emis- 
sion, the relative ratio of the integrated emission from up- 
stream and downstream is larger t han that calculated in the 
test p article case (see for instance IZirakashvili fc AharonianI 
toot) ) because particles upstream feel a smaller mean mag- 
netic field in the non linear case, hence the upstream syn- 
chrotron fiux is smaller than in the test particle case. The 
contribution of electrons to ICS from the upstream region 
is also smaller than the corresponding downstream contri- 
bution, although in the cutoff region the former provides an 
appreciable fraction of the total gamma ray fiux. 

We conclude the review of the results of our calculations 
with a discussion of the integrated spectrum of accelerated 
protons and electrons. These are shown (multiplied by p^) 
in Fig.|S]for both the hadronic (thick lines) and the leptonic 
scenario (thin lines). The break in the electron spectra is due 
to the onset of synchrotron losses in reaching equilibrium. 
In both scenarios the maximum energy of electrons is due 
to the balance between acceleration and synchrotron losses, 
although in the leptonic case the situation is borderline with 
the age of the system becoming comparable with the time 
for synchrotron losses. The most important comment on the 
spectra of accelerated particles, especially protons, is the 
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Figure 6. Spatially integrated spectra for accelerated protons 
(solid lines) and electrons (dashed lines) as functions of momen- 
tum. Thick and thin lines refer to the hadron dominated scenario 
(see Fig. [S} and lepton dominated scenario (see Fig. [JJ respec- 
tively. Circles represent the electron distribution as obtained by 
Tanaka et al. (2008). 



concavity of the spectra which appears to be barely visible. 
A standard non linear calculations with the shock parame- 
ters adopted here would lead to strongly modified (concave) 
spectra. The reduction of this effect is a crucial consequence 
of the dynamical reaction of the amplified magnetic field on 
the precursor, which leads to a reduction of the spectral con- 
cavity (see discussion in i|2.2|l . The cutoff in the spectrum of 
accelerated protons is due to finite acceleration time. Fig. [6] 
also provides a visual estimate of the value of Kep as it 
appears in Eq. \5\ namely Kep = 6 x 10~^ for the hadronic 
scenario and 1.4 x 10~^ for the leptonic one. Finally, in Fig. [6] 
we also compare our electron energy spectru m, as obtained 
in the hadronic scenario, with that derived bv lTanaka et al.l 
(2008) (see their Fig. 12) directly from the Suzaku X-ray 
data (circle points). The two spectra are in good agreement, 
with the possible exception of the few data points at the 
highest energies. 



4 CONCLUSIONS AND CRITICAL 
DISCUSSION 

The detection of high energy gamma rays from SNRs has 
long been considered as the most important tool to test the 
supernova paradigm for the origin of galactic cosmic rays. 
We now have a couple of positive gamma ray detections in 
the 1 — 100 TeV energy range, but it remains unclear whether 
this emission is of hadronic or leptonic origin. A hadronic 
origin would point to the acceleration of cosmic rays up to 
energies close to the knee, which is the main motivation for 
this work. 

We investigated the problem in the context of the non- 
linear theory of particle acceleration at SNR shocks, includ- 
ing the dynamical reaction of the accelerated particles, the 
generation of magnetic fields as due to streaming instability 
excited by cosmic rays, and the dynamical reaction of these 
self-generated fields on the plasma. We specialized our calcu- 
lations to the SNR RX J1713. 7-3946. Previous attempts to 
carry out similar calculations did not include the magnetic 
field amplification (which was input by hand so as to fit the 
data) and its dynamical reaction. The latter is crucial to re- 
duce the compression factors to the levels needed to explain 
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the few available observations (|Caprioli et al.|[2008al ). while 
in previous work on the topic, this proble m was faced by in- 
vokin g turbulent heating (see for instance (jBerezhko fc VoUg 

[iooi)). 

The non-linear calculations are applied to two cases 
with very different acceleration efficiency. In the first sce- 
nario, the shock accelerates a fraction ~ 10"* of the incom- 
ing particle flux, which results in large acceleration efficiency 
and effective magnetic field amplification through resonant 
streaming instability. This scenario provides a benchmark 
for a hadronic interpretation of the IfESS gamma ray obser- 
vations of RX J1713. 7-3946. In the second scenario a fraction 
~ 10"'' of the incoming particle flux through the shock gets 
accelerated. This serves as a model for a leptonic origin of 
the HESS data on RX J1713.7-3946. 

The main discrimination between the two scenarios re- 
sides in the strength of the magnetic field at the shock: effi- 
cient cosmic ray acceleration also leads to the generation of 
magnetic field upstream of the shock which by compression 
translates into a downstream field of B2 ~ 100/iG, compat- 
ible with the value which has been inferred from the inter- 
pretation of the morphology of the X-ray emission as due 
to synchrotron losses. The second scenario requires that the 
magnetic field is very weakly amplified (if at all), so that 
-B2 ~ 15 — 25/iG. In this case the X-ray morphology is in- 
terpreted as a result of damping of th e magnetic field in the 
downstream plasma l|Pohl et al.ll2005l ). 

In the following we discuss these two scenarios in terms 
of their problems and virtues, which also leads us to com- 
ment on the origin of the X-ray filaments. 

Our calculations show that for the parameters typically 
adopted for RX J1713. 7-3946, and assuming ^ = 3.7 for the 
injection, particle acceleration of protons works efficiently. 
X-ray and gamma ray data can be fit very well within this 
approach, and the strength of the downstream magnetic field 
turns out to describe well the brightness profile of the X- 
ray emission without further adjustments. The cutoff in the 
gamma ray emission leads to conclude that a cutoff in the 
parent proton spectrum is at ~ 10^ GeV. We stress that 
this is the cutoff in the proton spectrum which is being ac- 
celerated at the present time, while the highest energy which 
can be achieved in a SNR is reached at the very beginning 
of the Sedov phase. During the Sedov phase, particles which 
reached higher energies at previous times escape the system 
from upstream (and indeed these are the only particles which 
are allowed to leave the system from upstream). The shape 
of the X-ray spectrum is taken from the recent Suz aku ob- 
servations of RX J1713. 7-3946 (iTanaka et al.lbOOSi '). which 
lead to a cutoff energy in the X-ray spectrum requiring a 
shock velocity of 4300 km/s within the framework of our 
model. For Bohm diffusion the energy of the X-ray cutoff is 
very weakly dependent on the magnetic field strength in the 
shock region (in fact in linear theory it would have been in- 
dependent of B). In this scenario the ICS contribution from 
high energy accelerated electrons is negligible compared with 
the contribution of tt'' decays, unless an unreasonably large 
energy density of target IR light is adopted. 

The best fit is achieved by using a density of the up- 
stream gas of no = 0.12 cm~^. The implications of this re- 
sult are far reaching and deserve some discussion: as shown 
in Fig. |3] the emission expected from the gas behind the 
shock as a result of thermal bremsstrahlung emission of elec- 



trons exceeds the observed flux by a factor ~ 6. This conclu- 
sion cou ld be also derived based on simpler calculations, as 
done by 'Katz & Waxman' (2008!). However, as also argued 
by [Cassam-Chenai ct al . (2004), the electrons and protons 
in this SNR are not expected to be in thermal equilibrium 
because of the relatively low density of the m e dium . This 
possibility is also discussed bv lGhavamian iFaD(|2003), who 
find that for a remnant with velocity 3000 km/s the ratio 
of the electron to proton temperature is ~ 0.02 and scales 
as ~ V^^ (recall however that these conclusions were ob- 
tained in the context of linear calculations, and could be 
even stronger if the non-linear effects are taken into ac- 
count). We showed that such a temperature ratio would lead 
to predict a thermal X-ray emission well below the sensitiv- 
ity of Suzaku and ASCA. On the other hand, although the 
age of this remnant is too short to allow Coulomb scattering 
to equilibrate electrons and protons downstream, one should 
check that Coulomb scattering acting for ~ 1600 years does 
not heat electrons enough (about 1 keV) to excite lines in 
iron atoms, which would then become easily visible on top 
of the synchrotron emission (see for instance lEUison et al.l 
(2007) for a discussion of this effect). 

iKatz fc WaxmanI (120081 ) also criticized the hadronic sce- 
nario based on the lack of an excess g amma ray emission in 
the d irection in which, according to ICassam-Chenai' et al.l 
Hooi) the shock appears to have impacted a nearby molecu- 
lar cloud, where the target density for pp collisions should be 
higher. As simple as this inference may sound, it is not that 
straightforward: the gas in the molecular cloud is expected to 
be mostly neutral, a fact that has at least two important con- 
sequences. First, only the ionized component takes part in 
the acceleration process (recall that in order to be scattered 
back towards the shock from downstream, incoming parti- 
cles should be ionized within a few times the shock thick- 
ness). This implies that the efficiency of particle acceleration 
in the proximity of the molecular cloud could be much lower 
than in less dense parts of the surrounding medium. More- 
over, in a mostly neutral medium, the Alfvenic turbulence 
on which efficient particle acceleration relies can be effec- 
tively damped, thereby leading to a lower maximum energy 
of accelerated particles. These two instances show that hav- 
ing a dense target for pp collisions does not necessarily imply 
that there should be a larger gamma ray flux from tt" decay. 

As already stressed, the actual trademark of the 
hadronic scenario discussed above is the prediction that the 
downstream magnetic fleld in RX J1713. 7-3946 is B2 ~ 
100/iG, which would imply that the filaments observed 
in X-rays are d ue to sever e sync hrotron losses of accel- 
erated electrons. iPohl et all (|2(305l ) argued that a narrow 
region of X-ray emission could also be due to damping 
of the magnetic field down stream of the sho c k. Th e argu- 
ment was also adopted by iKatz fc WaxmanI (|2008l ) in or- 
der to suggest that HES S gamm a rays must be of leptonic 



origin. Katz fc WaxmanI (|2008l ) also quote the results of 



iRothenflug et aL I (|2004l ) on SN1006 to stress that filamen- 



tary structures are observed also in the radio and not only in 
X-rays, therefore the filaments cannot be d ue to synchrotron 
l osses . It should be noticed, however, that IRothenflug et al.l 
use their data to rea ch the opposite conclusion. In 
fact IRothenflug et al.l (|2004l 'l plot the spatial proflles of the 
radio and X-ray emission, which show that the radio emis- 
sion profile has only a shallow bump, when compared with 
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the pronounced peak in the X-ray brightness. We should re- 
call that the thickness o f the damping region, as correctly 
proposed by IPohl et all ()2005h . is determined by using the 
convection equation for the waves, with no generation terms 
and only a damping term. This trivially leads to a solu- 
tion which is an exponential in the spatial coordinate with 
a cutoff scale which defines the damping scale. It is worth 
reminding that the exponential shape is the very reason 
why one should expect a filament. In order to avoid that 
an e qually pronounced filamentary structure appears in ra- 
dio, iPohl et all l|2005l ) are forced to assume that the field 
has a residual large scale component which is not affected 
by the damping: the radio emission then comes from lower 
energy electrons propagating in this residual field. This solu- 
tion appears to us rather fine tuned, a limitation that should 
be kept in mind when discussing the damping explanation 
of the X-ray morphology. This limitation is probably even 
more serious in the case of RX J1713. 7-3946, where the X- 
ray emitting region is comparatively broader and it is not 
ve ry easy to explain it with the damping models proposed 
bv lPohl et al] (120051 ) (see Eqs. (7) and (12) of their paper), 
although there are many parameters in their formulae that 
can be appropriately tuned. 

In our opinion, a rather serious concern about a 
hadronic origin of the HESS gamma ray emission, and in- 
deed about this remnant being an efficient cosmic ray accel- 
erator, is the low value of K^p predicted by the non-linear 
theory. Our calculations lead to Kep ~ 10"*, which is much 
lower than the value inferred from direct measurement of the 
electron spectrum in the 10 GeV energy region, at the Earth. 
At first sight one should conclude that if the gamma radia- 
tion detected by HESS is of hadronic origin, then electrons 
should have a different origin, but a dedicated investigation 
of the problem should be planned to address this issue, since 
the very problem of electron escape from the remnant needs 
to be carefully addressed. One obvious warning against easy 
conclusions in this respect is that the value of K^p measured 
at the Earth is not the same as inferred for a single SNR at 
a given age. It is rather an average over the temporal evo- 
lution of a remnant and a sum over the many SNRs that 
contribute a flux at any given time. In fact Kep is in general 
a function of time for a SNR, and it is possible that low en- 
ergy electrons are more effectively produced at later times 
when proton acceleration and magnetic field amplification 
become less efficient, and probably there is no appreciable 
gamma ray emission due to pp scattering. 

Now we discuss our attempt to fit the HESS data with 
a leptonic model. We stress that the calculations carried out 
here are exactly the same as described above, with the only 
difference that we assume a larger value for the injection 
parameter ^, chosen here to be ^ = 4.1 (corresponding to 
a fraction of accelerated particles ~ 10~^). The generation 
of magnetic field by streaming instability is also kept in the 
calculations, although the predicted field strength is clearly 
much lower than in the previous case. 

As shown in Fig. |4l a fit to the combined HESS and 
Suzaku data is possible if the highest energy HESS data 
points are ignored. The cosmic ray acceleration efficiency 
corresponding to this solution is ~ 2%. The magnetic field 
found as a result of streaming instability and further com- 
pression in the precursor and at the subshock is ~ 20/iG 
and the cutoff energy of the electrons is still determined by 



energy losses. This is an important point, because the cut- 
off in the electron spectrum is super-exponential (namely 
oc exp[—{E / Ecut)'^]) only if the rate of losses is oc E'^. If the 
maximum energy of electrons is determined by the finite 
age of the accelerator, then an exponential cutoff should be 
expected. The two situations refiect in profoundly different 
X-ray spectra in the cutoff region. As seen in Figs. |3] and [l] 
the Suzaku data are well fit by an electron spectrum with a 
super-exponential cutoff. 

The fit illustrated in Fig. 3] has two problems, typical of 
leptonic models: 1) the nearly fiat region in the HESS data 
cannot be well fit by using C MB photons al o ne as the target 
for ICS, as also discussed by iTanaka et all (|2008l ). In order 
to solve this problem, a population of IR photons needs to be 
assumed, with energy density ~ 20 times larger than in the 
ISM. We consider this as a problematic aspect of the leptonic 
model, since this population of photons should be introduced 
ad hoc, without a proper independent physical motivation. 
2) The highest energy data points in the HESS data cannot 
be fit with ICS of electrons. These data points have been 
curiously ignored in some of the previous literature on the 
topic, probably because of their relatively large error bars. It 
is however not easy to imagine how the highest energy data 
point, at £ ~ 100 TeV may disappear. Future observations 
of this remnant in the multi-TeV range will tell us if this 
might be the case. 

Figs. 13] and |3] clearly show that a major discrimina- 
tion power between the two scenarios can be achieved with 
GLAST. In the region 1 — 100 GeV, the two scenarios imply 
a flux which is different by about one order of magnitude. 
A different conclusion was reached bv lYama zaki (2008), due 
to an oversimplified treatment of particle acceleration in the 
non-linear regime, which ignored several effects (magnetic 
dynamical reaction, turbulent heating) leading to much less 
concave spectra of accelerated protons. 
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